COMPLEXITY OF RANDOM SMOOTH FUNCTIONS ON COMPACT MANIFOLDS. II 



LIVIU I. NICOLAESCU 

ABSTRACT. To each even, nonnegative Schwartz function w(t) on the real axis we associate a random 
function on a compact Riemann manifold (M, g) of dimension m. We investigate the expected distri- 
bution of critical values of such a random function and the behavior of this distribution as we rescale 
w(t) to w(et) and then let e — > 0. We prove a central limit theorem describing what happens to the 
expected distribution of critical values when the dimension of the manifold is very large. Finally, we 
explain how to use the e — > behavior of the random function to recover the Riemannian geometry of 
(M,g). 
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Notations 

(i) For any set S we denote by \S\ € Z>o U {00} its cardinality. For any subset A of a set S we 
denote by I a its characteristic function 

/,:^ { o,i>, 

(ii) For any point x in a smooth manifold X we denote by 5 X the Dirac measure on X concen- 
trated at x. 

(iii) For any smooth manifold M we denote by Vect(M) the vector space of smooth vector fields 
on M. 

(iv) For any random variable £ we denote by E(£) and respectively var(^) its expectation and 
respectively its variance. 

(v) For any finite dimensional real vector space V we denote by V v its dual, V v := Hom( V, R). 

(vi) For any Euclidean space V we denote by Sym(V r ) the space of symmetric linear operators 
V — > V. When V is the Euclidean space R m we set Sym m := Sym(R m ). We denote by 
l m the identity map R' m ->■ R m . 

(vii) We denote by S(R m ) the space of Schwartz functions on R m . 

(viii) For v > we denote by y L , the centered Gaussian measure on R with variance v, 

1 

j v (x)dx = e \dx\. 

\/2irv 

Since lim^oT^ = 5o> we set 7o := <^o- F° r a rea l valued random variable X we write 
X G N(0, v) if the probability distribution of X is j v . 
(ix) If \i and v are two finite measures on a common space X, then the notation \i oc 1/ means 
that 

1 1 

— ; — rA* = — ; — 

1. Overview 

1.1. The setup. The present paper is a natural sequel to [24, 25]. Suppose that (M, g) is a smooth, 
compact, connected Riemann manifold of dimension m > 1. We denote by |dV^| the volume density 
on M induced by g. We assume that the metric is normalized so that 

vol s (M) = 1. (*) 

For any u, v G C°°(M) we denote by (u, v) g their L 2 inner product defined by the metric g. The 
L 2 -norm of a smooth function u is denoted by ||tt||. 

Let Ag : C°°(M) — > C°°(M) denote the scalar Laplacian defined by the metric g. Fix an or- 
thonormal Hilbert basis (^k)k>o of L 2 (M) consisting of eigenfunctions of A g , 

A g V k = \ k V k , ||* fc ||=l, k < h A feo < A fel . 

Fix an even Schwartz function w : M — > [0, 00). For e > we set 

u, e (i) : = Vi G R. 

Consider random functions on M of the form 

w e = ^-Ufe*fc, (1.1) 

fc>0 
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where the coefficients u k are independent Gaussian random variables with 

E(u k ) = 0, var(u k ) = w £ (yf\~ k ). (1.2) 

Note that 

A\ = ^Af«A, ViV>0. 

fc>0 

The fast decay of w, the Weyl asymptotic formula, [10, VI.4], 

_2_ 

Afc ~ const. km , k —> oo, 
and Chebyshev's inequality imply that 

E P (^ ^ 1737 ) < £*^«WAid < oo. 

k>l k k>l 

Invoking the Borel-Cantelli Lemma we deduce that for any N the function A N u £ is almost surely 
(a.s.) in L 2 . In particular, this shows that u £ is a.s. smooth. 

The covariance kernel of the Gaussian random function u £ is given by the function 

S £ : M x M R, 
£ E (p,q) = E(u E ( P )u e (q)) = ^MV^k)^k(p)^k(q)- 

k>0 

The eigenfunctions ^ k satisfy the known pointwise estimates (see [20, Thm. 17.5.3] or [26, Thm 
1.6.1]), 

( m + v \ 
X k 2 J as k -)• oo, > 0. 

Since w £ is rapidly decresing the above estimates imply that <f £ is a smooth function. More precisely, 
<§ e is the Schwartz kernel of the smoothing operator 

w(eV&) : C°°(M) -> C°°(M). 

The asymptotic estimates in Proposition 2.2 show that the random field du £ satisfies the assumption 
in [1, Cor. 11.2.2] for Invoking [1 , Lemma 11.2.11] we obtain the following technical result. 

Proposition 1.1. The random function u £ is a.s. Morse ife<^Ll. □ 

For any u G C 1 (M) we denote by Cr(tt) C M the set of critical points of u and by D(u) the set 
of critical values 1 of u. 

To a Morse function u on M we associate a Borel measure /jl u on M and a Borel measure a u on 
M defined by the equalities 

peCr(tt) teM 

Observe that 

supp^ u = Cr(ii), supp<r u = D{u). 

When u is note Morse, we set 

L^u '■= \dV g \, cr u = Sq = the Dirac measure on E concentrated at the origin. 



The set D(u) is sometime referred to as the discriminant set of u. 
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Observe that for any Morse function u and any Borel subset BcS the number a u (B) is equal to the 
number of critical values of u in B counted with multiplicity. We will refer to cr u as the variational 
complexity of u. 

To the random function u e we associate the random measure <j Ue . Its expectation 

a £ = E{a Ue ) 

is the measure on R uniquely determined by the equality 



/ f(t)a £ (dt)=E( [ f(t)d<r U£ (dt)) , 



for any continuous and bounded function / : R — > R. In §2.1 we show that the measure a e is well 
denned for e <C l.We will refer to it as the expected variational complexity of the random function 
u £ . The main goal of this paper is to understand the behavior of a £ as e — >■ 0. 

1.2. Statements of the main results. Observe that if u : M — > R is a fixed Morse function and c is 
a constant, then 

Cr(c + u) = Cr(u), fi c+u = fi u , 

but 

D(u + c) = c + D(u), cr u+c = 5 c *a u , 

where * denotes the convolution of two finite measures on R. More generally, if X is a scalar random 
variable with probability distribution u x , then the expected variational complexity of the random 
function X + u is the measure 

E(a x+U ) = v x * CTu- 
lf u itself is a random function, and X is independent of u, then the above equality can be rephrased 
as 

E(a x+U ) = vx* E(a u ). 
In particular, if the distribution v x is a Gaussian, then the measure E(cr u ) is uniquely determined by 
the measure E(a x+U ) since the convolution with a Gaussian is an injective operation. It turns out 
that it is easier to understand the statistics of the variational complexity of a perturbation of u £ with 
an independent Gaussian variable of cleverly chosen variance. 

To explain this perturbation we need to introduce several quantities that will play a crucial role 
throughout this paper. We define 



-r — / w(\x\)dx, 



(27T) 

1 

(2V) 
1 

(2V) 

If we set 

I k (w) := I w(r)r k dr, (1.4) 



hm := /o \m / x\xlw{\x\)dx. 



then we deduce from [ 5, Lemma 9.3. 10] 

dA(x) I m -i(w) = 
>\x\=i J Ml) 



(2ir) m s m = ( [ dA(x) ) I m -i{w) = ^L t^w), 
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J lxl= /i^dA(x) j i m+l ( w ) = 2r(2 + f) WM = m(m + 2)r(f) WM. 

We set 

m ' (i^ m + 2 I m+1 (w;) 2 
The Cauchy inequality implies that I m+1 (u?) 2 < / rn _i(w;)/ m+ 3(u;) so that 

Tfl 

q m > — • (1-5) 

m + 2 

The sequence (g m ) m >i can be inteipreted as a measure of the tail of w, the heavier the tail, the faster 
the growth of q m as m — > oo; see Section 3 for more details. The quantity s m has a simple spectral 
interpretation (B.l), more precisely 

f g £ (p,p)\dV g (p)\ =tvw(eVA) ~s m e- m , e \ 0. 

We can roughly interpret the above estimate as saying that the random function u £ belongs to a 
Gaussian-Hilbert space U e C C°°(M) of dimension s m e~ m . 
We set 

r n : = max(l,g n ), 

and define u m > via the equality 

r « = ^ • ( L6 ) 

Set s m := s m + u m so that 
Observe that 

cj m = 0<==^ m = r m > l^=>s m = s m , (1.8) 
while the inequality (1.5) implies that 

Um ^ = o, lim ^ = 1. (1.9) 

m— >oo s m rn— >oo q m 

Choose a scalar Gaussian random variable X w r e \ with mean and variance u>(e) := u) m e~ m inde- 
pendent of u £ and form the new random function 

We denote by & £ the expected variational complexity of ii e . We have the equality 

* e = 7u,( £ ) * ° E , oj(e):=u m e~ m , (1.10) 

Note that 



N e = [ & £ {dt) = [ a £ {dt) 



is the expected number of critical points of the random function u e . 

To formulate our main results we need to briefly recall some terminology from random matrix 
theory. 
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For v E (0, oo) and TV a positive integer we denote by GOE^, the space Symjy of real, symmetric 
N x N matrices A equipped with a Gaussian measure such that the entries a^- are independent, 
zero-mean, normal random variables with variances 

var(au) = 2v, var(aij) = v, VI < i < j < N. 

Let pn,v '■ R — > R be the normalized correlation function of GOE^. It is uniquely determined by 
the equality 



/ f(X)p N , v (\)d\ = ±E GOE « N (tTf(A)), 



for any continuous bounded function / : R — > R. The function pN,v(ty a l so has a probabilistic 
interpretation: for any Borel set B C R the expected number of eigenvalues in B of a random 
A e GOE^r is equal to 



N / p N)V {X)d\. 
Jb 

For any t > we denote by Jl t : R — > R the rescaling map R 9 x i— > tx G R. If /i is a Borel 
measure on R we denote by (3^)*// its pushforward via the rescaling map The celebrated Wigner 
semicircle theorem, [3, 22], states that as N — > oo the rescaled probability measures 

(01 i )♦( pj\r,„ (A)dA) 
converge weakly to the semicircle measure given by the density 

1 Jv^-A 2 , |A|<v^ 

We can now state the main results of this paper. 
Theorem 1.2. For v > arcc? iV G Z>o we se? 

#A^( X ) := PN,v(x)e ±1 ^. 

(a) There exists a constant C = C m (w) that depends only on the dimension m and the weight w such 
that 

N £ ~ C m (w)e- m (l + 0(e)) as e ^ 0. (1.11) 

More precisely 

m 

Cm (w) = 2^A (^|J 2 r (^p-j J^ rm -i * e+ +1 , r J(yhi(y)dy. (1.12) 



(b) As e \ f/je rescaled probability measures 

1 



converge weakly to a probability measure cr m on R uniquely determined by the proportionalities 

cr m oc (7r m -i *^+i, rm (y))7i(y)rfy (1.13a) 

3C 0~ ! _L_ * Tln^i (y)dl/- (1.13b) 
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Corollary 1.3. As e — > the rescaled probability measures 

converge weakly to a probability measure ct m uniquely determined by the convolution equation 

l^m, * <X m = <X m . 



Theorem 1.4. Suppose that r m has a limit r £ [1, oo] as m — > oo. Then 



lim a m = 7 r+i 

m— >0O r 



Corollary 1.5. As m — > oo we /jare 



6 m M — ^=r 



/nm \ 2 / \ ird m 



m + 3\ ( 2I m+3 (w) 



trm \ 2 J \ir(m + 2)/ m+ i(u>) 



')?( 



(1.14) 



Following [1 , §12.2] we define symmetric (0, 2)-tensor h £ on M given by 

m+2 

= — — E{Xu £ (p),Yu £ (p)), VpeM, X,YeVect(M), (1.15) 



where Xu denotes the derivative of the smooth function u along the vector field X. 

Theorem 1.6 (Approximation Theorem), (a) For s > sufficiently small the tensor h £ defines a 
Riemann metric on M. 

(b) For any vector fields X,Y on M the function h £ (X, Y) converges uniformly to g(X, Y ) as e \ 0. 

(c) The sectional curvatures on h £ converge to the corresponding sectional curvatures ofgase\0. 



The construction of the metrics h £ generalizes the construction in [ ]. Note that for any g > we 
have a smooth map H £ : M — > L 2 (M, g) 

/ £ m+2 \ | i 

MB P ^S £ (p) := [-T—J 2^^ £ (VAfc) 2 *ifc(p)*fe G L 2 (M, 5 ). 
^ m ' fc>0 

Then /i e is the pullback by 3 e of the Euclidean metric on L 2 (M, g). Let us point out that [6, Thm.5] 
is a special case of Theorem 1.6 corresponding to the weight w(t) = . 

Theorem 1.6 coupled with the results in [28] imply that the metrics h £ converge C 1,a to g as 
e \ 0. The convergence of sectional curvatures coupled with the technique of harmonic coordinates 
in [2, 28] can be used to bootstrap this convergence to a C°° convergence. 

We should add a few words about the refined analytic result hiding behind Theorem 1.6. Fix a 
point p G M and normal coordinates (x l ) at p. The techniques pioneered by L. Hormander [18] 
show that as e \ we have the 1-term asymptotic expansions 

E(d 2 x ^u £ (p) ■ d 2 x3xJ u £ {p)) =h m e-( m+4 ^l + 0(e 2 )), (1.16a) 
/';(<•';„. " (p) • dl xJ u £ {p))=h m e-^( K l + 0{e 2 )). (1.16b) 



s 
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These and several other similar 1-term asymptotic expansions involving the Schwartz kernel of w e (VK) 
(see Proposition 2.2) are responsible for Theorem 1.2. All these 1-term expansions are independent 
of the background metric g. Note that (1.16a) and (1.16b) imply the estimate 

"{<>:,■<> (p) ■ i>;., « (p)) !<(;>:..,,■<> (p) ■ d 2 xixj u% P )) =o( £ -(-+ 2 )). (i.i?) 

Theorem 1.6 is equivalent with the following estimate 

E[dl ixi u%p) ■ dl jxj u%p))-E(d? ixj u%p) ■ d 2 xixj u%p)) ~d m ^.(p) £ -(-+ 2 ), 

where K'f- (p) denotes the sectional curvature of g at p along the 2-plane spanned by d x i , d xJ which 
is an obvious refinement of (1.17). In fact the proof of Theorem 1.6 is based on refinements of the 
1-term expansions (1.16a) and (1.16b) to 2-term expansions. 

The Schwartz kernel of w e {\J~K) has a complete asymptotic expansion as e \ (see [29, Chap. 
XII]) and Theorem 1.6 shows that the metric g becomes visible, and it is completely detected by the 
second order terms of this expansion. 

The convergence of the metrics h e has a pleasant consequence. Assume that M is oriented and 
m = dim M is even. To a Morse function / we associate the signed measure 

"/= E (-i) ind(/ ' p) v 

d/(p)=0 

where ind(/, p) denotes the Morse index of the critical point of the Morse function /. It is convenient 
to think of Vf as a degree 0-current. The random function u e determines a random 0-current v u e. 
It turns out (see Section 4) that the expectation of this current is a current represented by a rather 
canonical top degree form, namely the Euler form e^e (M) defined by the metric h £ which appears in 
the Gauss-Bonnet theorem. The Morse inequalities imply 

/ e h e(M) = f v e = E ( [ u u s) = x(M), 
Jm Jm \Jm ) 

and as a bonus we the Gauss-Bonnet theorem for the metric h £ . Letting e — > we obtain the Gauss- 
Bonnet theorem for g since h £ — » g and e^e(M) — > e g (M). In particular, this shows that v e 
converges in the sense of currents to e g {M), the Euler form determined by the metric g. 

1.3. A few comments. The remainder of the paper is organized as follows. Section 2 contains the 
proofs of the main results. In Section 3 we discuss various concrete weights w. In particular, these 
examples show that the limit r = limm^oo r m that appears in the statement of Theorem 1.4 can have 
any value in [l,oo]. Section 4 contains the details of the probabilistic proof of the Gauss-Bonnet 
theorem outlined above. 

To smooth the flow of the presentation we gathered in Appendices various technical results used in 
the proofs of the mains results. In Appendix A we describe the jets of order < 4 along the diagonal of 
the square of the distance function dist 9 : M x M — > R. These play a role in the spectral estimates 
of Appendix B and we could not find any reference describing these jets with the level of specificity 
we need. 

The technique of obtaining the spectral estimates in Appendix B is well known to experts in the 
field, but again, we could not find any reference describing these estimates with the level of precision 
required by our concrete investigation. In particular, we believe that Proposition B.5 is new. This 
result is an explicit two-step asymptotic expansion of the Schwartz kernel of w e (v~K) as e \ 0. 

In Appendix C we describe a few facts about Gaussian measures in a coordinate free form suitable 
for our geometric purposes. Finally, in Appendix D we have collected some facts about a family of 
Gaussian random symmetric matrices that appear in our investigation. 
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2. Proofs 

2.1. A Kac-Rice type formula. The key result behind Theorem 1.2 is a Kac-Rice type result which 
we intend to discuss in some detail in this section. This result gives an explicit, yet quite complicated 
description of the measure a 6 . More precisely, for any Borel subset B C R, the Kac-Rice formula 
provides an integral representation of cr e (B) of the form 

<r £ (B) = [ f £>B (p)\dV g (p)\, 

JM 

for some integrable function # : M — > R. The core of the Kac-Rice formula is an explicit 
probabilistic description of the density f £> B- 

Fix a point p G M. This determines three Gaussian random variables 

u £ (p) G R, du £ (p) G T;M, Hess p (-u e ) G Sym(T p M), (RV) 

where Hess p (it £ ) : T p M x T p M — > R is the Hessian of m u at p defined in terms of the Levi-Civita 
connection of g and then identified with a symmetric endomoiphism of T p M using again the metric 
g. More concretely, if (x l )i<i< m are g-normal coordinates at p, then 

m 

Hessp(w £ )^ = ^d 2 xlx3 u £ {p)d xl . 

i=l 

For e > sufficiently small the covariance form of the Gaussian random vector dii £ (p) is positive 
definite; see (2.3). We can identify it with a symmetric, positive definite linear operator 

S(du £ {p)) : T p M -> T p M. 

More concretely, if (x*)i<j< m are ^-normal coordinates at p, then we identify S(dii £ (p) ) with a 
m x m real symmetric matrix whose (i, j) -entry is given by 

Sij(du £ (p)) = E(d Xi u £ (p) ■ d xJ u £ (p)). 

Theorem 2.1. Fix a Borel subset Bci For any p G M define 

/ e , B (p) := (det(2vr5(« e (p) | det Hess p (« e )| • J B (w £ (p)) | ^ e (p) = o), 

where E ( var | cons ) stands for the conditional expectation of the variable var given the con- 
straint cons. Then 

& £ {B)= [ f e , B (p)\dV g (p)\. (2.1) 

JM 

a 

This theorem is a special case of a general result of Adler-Taylor, [1, Cor. 11.2.2]. Proposition 2.2 
below shows that the technical assumptions in [1, Cor. 11.2.2] are satisfied if e <C 1- 

For the above theorem to be of any use we need to have some concrete information about the 
Gaussian random variables (RV). All the relevant statistical invariants of these variables can be 
extracted from the covariance kernel of the random function ii £ . 
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2.2. Proof of Theorem 1.2. Fix e > 0. For any p S M, we have the centered Gaussian random 
vector 

( u £ (p), du £ (p),Ress p (u £ ) ) 6l(B T*M © Sym(T p M). 

We fix normal coordinates (x l )\<i< m at p and we can identify the above Gaussian vector with the 
centered Gaussian vector 

(u £ (p), (d x iU £ (p) )i<i< m , d 2 xix3 {u £ {p) )i<i,j< m ) el©« m © Sym m . 

The next result is the key reason the Kac-Rice formula can be applied succesfully to the problem at 
hand. 

Proposition 2.2. For any 1 < i,j,k,£ < mwe have the uniform in p asymptotic estimates as e \ 



E(u £ (p) 2 ) =s m£ - m (l + 0(e 2 )), (2.2a) 

E(d x ,u £ (p)d x3 u £ ( P ) ) = d m e-( m+ %-(l + 0(e 2 ) ), (2.2b) 

E(d 2 lx3 u £ (p)d 2 xkxe u £ (p) ) = h m e^ m+ ^ (5 tJ 5 M + SikSjt + 6 u 8 jk ) ( 1 + 0{e 2 ) ) , (2.2c) 

E{ u £ {p)d 2 xlx3 u £ (p) ) = -dme-^+^Sij ( 1 + 0(e 2 ) ) , (2.2d) 

E(u £ (p)d x ,u £ (p)) = 0(e- m ), E{d x ,u £ {p)d 2 xJx3 u £ {p)) = 0(e~^), (2.2e) 

w/zere s m = s m + u> m a«J f/ze constants s m , d m , h m are defined by (1.3). □ 



Proof. Denote by S £ the covariance kernel of the random function 

« £ = ^w(e) + u s- 

Note that 

<# e (p, g) = u(e) + S £ (p, q) = u m e~ m + £ £ (p, q). 
Fix a point p G M and normal coordinates at p defined in an open neighborhood Oo of p . The 
restriction of S £ to Oo x Oo can be viewed as a function £ £ (x, y) defined in an open neighborhood 

of (0, 0) in R m x R m . For any a, /3 € (Z> )| m we have 

E(d%u £ (p )d£(u e )) = d^d^ £ (x,y) x=y=0 . 
Proposition 2.2 is now a consequence of the spectral estimates (B. 1) in Appendix B. □ 

From the estimate (2.2b) we deduce that 

S( du £ (p) ) = d m e^ m+ ^ ( l m + 0(e 2 ) ) , (2.3) 

so that 

, m m(m+2) , . 

^{det S(u £ (p))\ = {d m )-e —{l + 0(e 2 )) as e -> 0. (2.4) 

Consider the rescaled random vector 

(s £ ,v £ ,H £ ) =:= (eTu^^H^diieip), e^ V 2 ii £ (p) ) . 

Form Proposition 2.2 we deduce the following (uniform in p) estimates as e \ 0. 

E((s £ ) 2 ) =s m (l + 0(e 2 )), (2.5a) 

E ( v £ v £ ) = dmSij ( 1 + 0(e 2 ) ) , (2.5b) 

E{H £ jH £ kl ) = h m (5 i3 5 M + <5 ifc <5 j£ + S u S jk ) ( 1 + 0(e 2 ) ) , (2.5c) 

B( s £ R £ 3 ) = -d m <% ( 1 + 0{e 2 ) ) , (2.5d) 

E(s £ vf) = O(e), E(v £ H £ k ) = 0(e). (2.5e) 
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The probability distribution of the variable s £ is 

1 x 2 
dls m {e){x) = 7= = . e 2 ^)\dx\, 
■ v /2vrs m (e) 

where 

s m (e) = s m + 0(e). 

Fix a Borel set Bel. We have 

£(|detV 2 u e (p)|I B (u e (p)) | du £ (p) = 0) = e — tf( | det # £ |I £ ™ B (s e ) | W £ = 0) 



m(m+4) C , P 2s m (e) 

■/ £(|detff £ | | s £ = x, v £ = 0) -. :\dx\ 

Je^B v /27TS m (e) 



(2.6) 



=: % , p ( e lTB) 

Using (2.4) and (2.6) we deduce from Theorem 2.1 that 

* £ {B) = e~ m 2 J M q E , p (efB)p L (p)\dV g (p)\, 

where p e : M — >■ R is a function that satisfies the uniform in p estimate 

p e (p) = 1 + O(e) as e -»• 0. (2.7) 

Hence 

m 

«- = (j^.) * »,(B) A (P)I*4(P)|. (2.8) 

To continue the computation we need to investigate the behavior of q £jP (B) as e. More concretely, 
we need to elucidate the nature of the Gaussian vector 

(H £ \s £ = x, v £ = 0). 

We will achieve this via the regression formula (C.3). For simplicity we set 

Y £ := {s £ ,v £ ) el®l m . 

The components of Y £ are 

Y £ = s £ , Y £ = v £ , \<i<m. 
Using (2.5a), (2.5b) and (2.5e) we deduce that for any 1 < i, j < m we have 

E{Y £ Y £ ) = S m Soi + 0(e), E(Y £ Y £ ) = ci m <% + 0(e 2 ). 

If S(Y £ ) denotes the covariance operator of Y, then we deduce that 

S(Y%} = ^-<S« + 0(s), SiY^- 1 = -LB* + 0{e). (2.9) 

We now need to compute the covariance operator Cov(H £ ,Y £ ). To do so we equip Sym m with the 
inner product 

(A,B) = tr(AB), A,B eSym m 
The space Sym m has a canonical orthonormal basis 

Eij, 1 < i < j < m, 
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dij 



where 

E =S Eij ' i = j 

and Eij denotes the symmetric matrix nonzero entries only at locations and and these 
entries are equal to 1. Thus a matrix A £ Sym m can be written as 

where 

v V2oij, i < j- 
The covariance operator Cov(H £ ,Ye) is a linear map 

Cov(H £ ,Y £ ) : R © R m —?■ Sym m 

given by 

Cov(H £ ,Y £ ) (f> a ej = ]T E(H £ 3 Y £ )y a E tj = £ E^Y^E^, 

\a=0 / i<j,ct i<j,ot 

where e , e l5 . . . , e m denotes the canonical orthonormal basis in R © R m . Using (2.5d) and (2.5e) 
we deduce that 

Cov(H £ ,Y £ ) lj2y a e a j = -y ^lm + O(e). (2.10) 

\q=0 / 

We deduce that the transpose Cov(H £ ,Y £ ) V satisfies 

Cov{H £ ,Y £ ) v (v^aijEi^j =-d m tr{A)e + O{e). (2.11) 

Set 

Z £ := (H £ \s £ = x,v £ = 0) - E(H £ \s £ = x,v £ = 0). 
Above, Z £ is a centered Gaussian random matrix with covariance operator 

S{Z £ ) = S{H £ ) - Cov{H £ ^SiY 6 )- 1 Cov(H £ ,Y £ f . 

This means that 

E{zl j r M ) = (E ij ,S(Z E )E k£ ). 
Using (2.9), (2. 10) and (2. 1 1) we deduce that 

Cov(H £ , Y £ )S(Y £ )~ 1 Cov(H £ , Y £ ) v [ fiy.Ey | = ^ tr(A)l m + 0(e) 



m 



E{(z £ 3 ) 2 )=h m + 0(e), E(z £ t z £ j ) = h m -^ + 0(e), Vi < j, 
£((4) 2 ) =3^-^ + 0(6), Vi 

and 



E(z £ jZ £ ke ) = 0(e), Vi<j, k<£, (i,j)^(k,. 
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We can rewrite these equalities in the compact form 

dl 



Note that 



E (z!j4e) = [h m - — ) Sij5 ki + h m (5 ik 5 je + 5 ie 5 jk ) + O(e). 



h _ ZUL 

1 "m. 



d 2 m (1.7) r m - 1 



We set 



so that 



(r m - 1) 

m. 



E i z !j z ke) = 2K m h m 5ij5 ke + h m (5 ik 5 je + S ie 5 jk ) + 0(e). 
Using (C.4) we deduce that 

E{H £ \s £ = x, v e = 0) = Cov{H £ , Y^SiY'y^xeo) = - ^l m + 0(e 



(2.12) 



We deduce that the Gaussian random matrix (H e \s £ = x,v £ = 0) converges uniformly in p as e — > 



to the random matrix ^4 - 
described in Appendix D. Thus 



l m , where A belongs to the Gaussian ensemble Sym„ 



m 



lim q £)P (B) = qoo(B) := / -Eg ym 2« m fe m ,fc m 



£->0 



B 



det(A-^l m )|) 



7TSr 



(/lm) * / ^SvmS,"- 1 ( I det ( A ~ x ~7T = '^ m ) I ) 



e 2 » 



B 



Sym 



V2 



-.dx 



TTS r 



y 

e 2 



where 



This proves that 



1 _ E Sym 2K, n ,i (\det(A- a m yl m ) | ) —f=dx 



2tt 



dm (1.7) 1 



y f s m h r 



l™\ m) -^ £ , P (S) = (^m)2 / £ g>> l 



B 



det ^ 



:1, 



_ a_ 
e 2 



dy . 



= :/j, m (B) 

Using the last equality, the normalization assumption (*) and the estimate (2.7) in (2.8) we conclude 



(X i ) & £ {B) = e 



-m I j ""rn 

2ird v 



H m (B) + 0{e) ase^O. (2.13) 



In particular 



h m 



2Tvd r , 



/i m (M) + 0(e) as e -»• 0. 



Observe that the density of /z m is 



s y m 



det A 



e 2 



2tt 



(2.14) 



(2.15) 
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(A = sfr^A) 



{2k m T m — T m 1) 

(D.7b) 



Sym,, 



r m 2 22 (2r m ) — r 

„m±± i /m + 3 
= 2 2 r£J 



detf A - yl T 



e 2 



2tt 



i / m + 3 



(7r m -l * )(y)7i(y)- 



(7r m -l * 6 'm+l,r m )(2/)7l(y)- 

This proves part (a) and (1.13a) in Theorem 1.2. To prove (1.13b) we distinguish two cases. 
Case 1. r„, > 1. From Lemma D.2 we deduce that 



-£/ri 2^771 ,1 

Sym m m ' 



det A 



:1, 



m+3 _ f m + 3 \ 1 

2— r ' 



where 



Pm+i,i(A)e ^ (m+ W + dA, 



_2 1 



(2.16) 



t : = 



Kvyi 1 r m + 1 



Thus 



d^m m+3 / m + 3 \ 1 (^+l-2r m )y 



dy 



2— r 



7TK r 



VCA-fr 2 +1) » ) 2 



m+3 ^ /m + 3\ 1 

2 2 r 



2 / y/2lTK r 

An elementary computation yields a pleasant surprise 



p m+ n(A)e ^ l ™ ^^dA. 



1 



4r 2 
Now set 



1 

A - y* 



2(r ro + l) 4 IV 2 ^™- 1 ) V2(r m -1) 



We deduce 



G^m „IZI+3 /m + 3\ 1 

= 2 2 r 1 



2 y 2tt,/k: 



(r m - 1)' 

p m+ iAX)e-^e-^ X -^ 2 d\. 



(A := V^A) 



m+3 _ / 771 + 3\ 1 

2~ r 



(73.6) m+3 /m + 3 

= Z 2 



\J K'rnTraf^n 



, , , , £m \2 r m 0m r\ ,,\2 

Pm+i,i/r m (A)e^ I g LA2 (i7 ^ i (y - A)dA. 



2 m±4 r /m + 3 



Pm+i,i/r ro (A)e r ™ A 'd 7 i (y - A)dA 

Pm r m 
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= 2 m i i T[ r ^— \ I p m+ i A / rm {\)e r ™ ^d^rm-i (y - X)d\ 
V z J Ju Tm 
Using the last equality in (2.13) we obtain the case r m > 1 (1.13b) of Theorem 1.2. 

Case 2. r m = 1. The proof of Theorem 1.2 in this case follows a similar pattern. Note first that in 
this = so invoking Lemma D.l we obtain the following counterpart of (2.16) 



E 



GOEi 



det( A - yl 7 



+4 f m + 3 



2— r 



e 4 Pm+i,i(y)- 



Using this in (2.15) we deduce 

dfl r , 



m+4 _ /m + 3 



2— r 



e 4 p m+1>1 (?/) 



which is (1.13b) in the case r m = 1. This completes the proof of Theorem 1.2. 
2.3. Proof of Corollary 1.3. According to (1.10) we have 



7u 



Thus 



Hence 



7 "m * ( 3? 1 I 



31 i 



1 



Km -Ji£rn * I 3? 1 



(T = rr, 



□ 



We can now conclude by invoking Levy's continuity theorem [21, Thm. 15.23(h)] or [30, Thm. 
2.4]. □ 



2.4. Proof of Theorem 1.4. We have 



1 _j_ * J rm-l dy 



(2.17) 



where 



J- 1 _ ' T-m 



e 4 , 



and 



K„ 



_ x j_ * 7 r m -i (v)dy 



± _i_ (A)dA = / p m+li ^(A)e 4 dA. 



We set 



^m(A) := p m+1) i(A), Roo{x) = ^i"{| x |<2}V / 4-x 2 . 
Fix c E (0, 2). In [24, §4.2] we proved that 



lim sup \Rm(x) - Roa(x)\ = 0, 

m ^°°|x|<c 

sup \R m {x) — Roo{x)\ = 0(1) as m — > oo. 

\x\>c 



(2.18a) 
(2.18b) 



Then 



/> m+ l,^i_(A) = \I^Rr, 





A ) , 0" ! (A) 





A I e 4 

m 



We now distinguish two cases. 
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Case 1. r = lim m _ 5 . 00 r m < oo. In particular, r G [1, oo). In this case we have 

K m = M.f R m (J^\) e-^dX, 
V m Jr \\ m J 

and using (2. 18a)-(2.18b) we deduce 

lim / R m (a/^X) e-^dX = R^O) [ e'^dr = R^O) J— . 



Hence 



Now observe that 



K m ~K'= Roo(0)\ — as m -> oo. (2.19) 
V m 



J-0- 1 , (A)dA = — l —-R m ( J^X) -^Le^dX 



;R m (\^X) 1 ^(dX) 



Roo (0) \\ m J r m 
Using (2.18a) and (2.18b) we conclude that the sequence of measures 

-^—Q- , ! (A)dA 

converges weakly to the Gaussian measure 72 . Using this and the asymptotic equality (2.19) in (2.17) 

7' 

we deduce 

lim & m = 72 * 7 r-i = 71+1 . 

171 — >00 r r r 

This proves Theorem 1.4 in the case r < 00 since 

w m (1.9) 

7 -> * <r m = cr m and iim — — = U. 
Case 2. linim^oo r m = 00. In this case we write 

! WdX = \[^R m L[^x) 1 ^{X)dX. 

Lemma 2.3. 77je sequence of measures 

\ y m / ,m 
converges weakly to the measure Roo(0)5q. 

Proof. Fix a bounded continuous function / : R — >■ R. Observe first that 

lim / f ,p^) - i?oo /(A) 7 j_(A)dA = 0. (2.20) 



=D m 

Indeed, we have 



\\\<c- 



--■DL 
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1 1 v" r «l 

S v ' 

= :D" 

m 

Observe that 

D' m < sup \R m (x) - Roo^l [ _ /(A)7j_(A)dA 

|x|<c ■ / I A I< C ^ 

and invoking (2.18a) we deduce 

lim D = 0. 

m— too 

Using (2.18b) we deduce that there exists a constant 5 > such that 

D' m <S [ 7^(A)dA. 

On the other hand, Chebyshev's inequality shows that 

2 



|A|>c 



Hence 



7 ^(A)dA< 



lim LC = 0. 

m— »oo 



This proves (2.20). 

The sequence of measures ^j2_ (\)d\ converges to 5q so that 

Roo(0)f(0)= lim I Roo(0)f(Xh^(X)d\. 
Using (2.20) and the above equality we deduce that the conclusion of the lemma is equivalent to 

lim / (r^W-R^ /(A) 7 _2_(A)dA = 0. (2.21) 

" v ' 

=Fm 

To prove this we decompose F m as follows. 



F m = ! ^_ -Roo(O) - -Roo \ X /(A) 7 _3_(A)dA 



|A| <m -i^M V V V m 



r. 




rm 



[ , ^ (Roo(0)-Roo( x ^x)) f(Xh^(X)dX. 



:F" 



Observe that 



F' m < sup \Roo(0) - Roo(x)\ [ 1 ^ /(A) 7 ^(A)dA 



and since i?„o is continuous at we deduce 



lim F' m = 0. 

m— >oo 
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Since and / are bounded we deduce that there exists a constant S > such that 

J\\\ >m -?4sL r ™ 



'|A|> 

On the other hand, Chebyshev's inequality shows that 



/ x _ 7^(A)dA< 



' ' s/rrn 

Hence 

lim F" = 0. 

m— >oo 

This proves (2.21) and the lemma. 
Lemma 2.3 shows that 

I Air 

K m ~ K' m = J—R oo (0), 
V m 



and 



On the other hand 



so that 



lim — 6~ 1 ! (A)dA = 5 . 



lim 7 m,-i (A)rfA = 7i(A)dA, 



lim <r m = 5 * 71 = 71. 

m— >-oo 

This completes the proof of Theorem 1 .4. 

2.5. Proof of Corollary 1.5. Using (2. 14) we deduce 

m 

m 

= ^rM fe- T /V ^(A) d A + 0(E ) 



/ \^' iu "m/ JR ' r m 



=Cm(w) 

Lemma 2.3 implies that as m — > oo we have 



Air 2 
#m~\/— -Roo(O) 



m \Jixm 
We deduce that 



m+6 , \ / 7 

2 2 _ / m + 3 \ / Or 
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2.6. Proof of Theorem 1.6. Fix a point p G M and normal coordinates (x l ) near p. The equality 
(2.2b) shows that as e — > we have the following estimate, uniform in p. 

E(d x iU £ (p)d x ms(p) ) = dm.£- (m+2) (Sij + 0(e) 2 ) ) . 

Hence 

h B (d xi , ) = Sij + 0(e 2 ) = g p (d xi , 8 xJ ) + 0(e 2 ). (2.22) 

This proves (a) and (b) of Theorem 1.2. 
With p and (x l ) as above we set 

d a+b S £ {x,y) 
dx %1 ■ ■ ■ dx la dy^ ■ • • dyi 



e>e ._ 

il,...,i a ;ji,...,j b ■ p,^ . . . Ftrr-iaFl^h . . . f) v j b \ X =V=0> 



hfj ■= h £ p (d xl ,d xJ ), 1 < i,j < m. 

We denote by Kf- the sectional curvature of h e along the plane spanned by d x i , d xl . Using [ , Lemma 
12.2.1] and that the sectional curvatures of a metric are inverse proportional to the metric we deduce 
as in [24, §3.3] that 

a .. g £ . . . 



V £ m+2 gege _Ige\2- 

Using Proposition B.5 we deduce that there exists a universal constant Z rn that depends only on m 
and w such that 

" ^;ij = e- (m+2) )Z m ^(p)( 1 + 0(e 2 ) ), 
where Kij(p) denotes the sectional curvature of g at p. The estimate (2.2b) implies that 

S ¥h ~ = d 2 m e- 2 ^ ( 1 + 0(e 2 ) ) . 

Thus 



Kt 3 = ^ Kl] {p)(l + 0{e 2 )). 



in 



To determine the constant ^ it suffices to compute it on a special manifold. Assume that M is the 
unit sphere S m equipped with the round metric. This is is a homogeneous space equipped with an 
invariant metric g with positive sectional curvatures. The metrics h e are also invariant so there exists a 
constant C e > such that h £ = C e g. The estimate (2.22) implies that C e = 1 and thus Kf- = Kij(p) 

so that % ta - = l. □ 

dm 

3. Some examples 
We want to discuss several examples of weights w. Observe first that 

/ \ D i x Im-l{w)I m+ z(w) 

r n (w) ~ R m (w) = as m ->■ oo. 

Moreover 

R n {w £ ) = R n (w). 

Example 3.1. Suppose that w(t) = e~ t2 . In this case <r> £ is the Schwartz kernel of the heat operator 
e~ £ A whose asymptotics as e — > have been thoroughly investigated. The momenta (1.4) are 

I k {w) = r t k e- t2 dt = X - f°° s k -^e~ s ds = l -T ( 
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Hence 

r(f)r(f + 2) m + 4 m(m + 4) 

so that r rn = 1 for all m. Moreover, in this case we have 



I m+ i(w) 

so that 



m + 2, 



m + 6 , 

2— /m + 3, 

^m(w) — m+1 T — — ) asm^oo, 

/mir 2 V z 



and Stirling's formula implies 



TfL 

log C m (w) ~ — log m as m — )• oo. (3.1) 

□ 



Example 3.2. Suppose that 

w(t) = exp(-(logt)log(logt)), Vi > 1. 

Observe that 

rl rco 

Ik(u>) = / r k w(r)dr + / r fc exp(— (logr) log(log r) )dr. 

This proves that 

/oo 
r fc exp(— (log r) log(log r) )<ir as /c — > oo. 

Using the substitution r = e* we deduce 

/>oo 

J fc = / e ( fc+1 )'-' 10 ^dt. 

J 

We want to investigate the large A asymptotics of the integral 



;>oo 

T x = e~ Mt) dt, <f> x (t) = Xt - tlogt. (3.2) 
J o 

on the Laplace method [9, Cha 

&(i) = A-logt-l, 



We will achieve this by relying on the Laplace method [9, Chap. 4]. Note that 

1 

T 

Thus 4>\{t) has a unique critical point 

T = r(A) —e^ 1 . 
We make the change in variables t = ts in (3.2). Observe that 

Ae A_1 s - e A_1 s log^- 1 ^) = e A_1 s - (A - l)e A_1 s - e A_1 log s = e A-1 s(l - log s 
and we deduce 



poo 

T x = t e - Th{s) ds, h(s) = s(logs-l). 



The asymptotics of the last integral can be determined using the Laplace method and we have, [9, 
§4.1] ' 
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Hence 

J k = T k+1 ~ y/2iTT(k + l)e T(fc+1) = V2^e ek as k -»• oo. 

In this case 

R m {w) — > oo as m — >■ oo. 

Note that 

fom _ 2I m+3 (w) 

d m (m + 2)/ m+ i O) ' 

We deduce that 

g — ~ e — e as to — > oo. 



Hence 



logC m (u;) ~ ye m+2 (e 2 - 1) as m -> oo. 



Example 3.3. Suppose that 

w(r) = exp(-C(logr) a ), C > 0, r > 1, a > 1. 
Arguing as in Example 3.2 we deduce that as — >■ oo 

/OO /'OO 
r fc exp( -C(logr) a )dr = J e {k+1)t - cta dt. 

Again, set 

/•oo 

T\ := / e'^dt, cp x (t) := Ct a - At. 
J o 

We determine the asymptotics of T\ as A — > oo using the Laplace method. Note that 

(f>' x (t) = aCt a ~ l - A. 
The function <fi\ has a unique critical point 



Observe that 



aC 



( A \~ _j_ 
b x (Ts) = a{s a -bs), a := - - , 6:=a»-i. 



/■oo 

T A = r(A) / e- a ( sQ - fes )ds. 

JO 



We set g(s) := s a — bs. Using the Laplace method [9, §4.2] we deduce 



T\ ~ r(A)e"^ 1 ) = ./ 2n e «ft-D. 

V o/(l) V aa(a - 1) 

Hence 

/A Q \ - 1 

logT A ~ ( — ^ =: z(a,C)A— . 

Hence 

log R m (w) ~ log T m + log T m+4 - 2 log T m+2 
~ Z(a, C) ( m^r + (m + 4)^r - 2(m + 2)^r 
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= Z(a,C)m^ ( 1+ (l + -) a - 1 - 2(l + - 
\ \ mJ \ m 

_s_ 8 a ( a \ 8aZ(a) 2=a 

~ Z(a, C)m-i x — x - 1 = . 

m z a — 1 \a — 1 / (a — l) z 

Hence 

{00, a < 2, 

e 16Z(2,C) a = 2 
1, a>2. 

which shows that r can have any value in [1, 00]. Note that in this case 

log I m+3 (w) - log I m+ iH ~ Z(a, C)m^ ( (l + -) ~' - (l + - ^ 

2Z(a,C) _JL_ 

~ m^- 1 , m — > 00, 

a — 1 

so that 

. „ , s Z{a,C) 
a — 1 



□ 



Example 3.4. Suppose now that w is a weight with compact support disjoint from the origin. For 
example, assume that on the positive semi-axis it is given by 



I e \x-c\ < 1 -1 

7i>(XJ = < 1 , C > 1. 

I 0, \x - c| > 1, 



Then 



c+l 1 rl 



I k (w)= I t k e i-c-^dt = / (t + c) fc e T^di 



it 



Observe that 
On the other hand 

and we deduce 



lim c~ k lZ = 0. 
Jo 

c k [ e~^dt < /+ < (c+l) k [ e~^dt. 
Jo Jo 

Hence the asymptotic behavior of Ik{w) is determined by it. We will determine the asymptotic 

behavior of it by relying again on the Laplace method. Set a := (c + 1) so that 

rl 1 r— 1 poo 2 

1+ = / ( a - tfe'^dt = a k / " (1 - s) fc e"^ds = a k (u - l) k u~ {k+2) e'^ 1 du. 

Jo Jo Ja 

Consider the phase 

0fi(s) = - bg(l - S) 5-3 , fr\0, 
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and set 

P h = a 1 h f° e ^W 
Jo 

so that 

It = Pi/k- 

We have 

The phase as a unique critical point r = r(7t) G (0, l/o) satisfying 

a 2 r 3 a 2 r 3 



so that 



Set 



n = 2(1-0 + 



T = I ^ j 3 (l + 0(hk ) ) as h \ 0. (3.3) 



1 a 2 r 4 2ft 3 1 /2fr 2 y 
v:=v(h ) : =-—^~ _^^J_= _ . (3.4) 



t 1/a — r 



We make the change in variables s = r + y^x and we deduce 

JJ{h) 
We claim that 

lim / e <t>h{T+^x)-4, h {T) dx = (' e -4dx = \P2^. (3.5) 
>>-> Jj(h) Jr 
It is convenient to think of r as the small parameter and then redefine 

o 2 t 3 
^ 2(1 -r) 
and think of v as a function of r. Finally set o := yjv and 

ip T (x) := <f) h ( T ){T + ox) - ft ( r )(r) = 2 ^2 T 3 T ^ 1 og( 1 - s) - ^2 

2(1 -r)/ , . , x \ 1 / 1 1 

' log(l - r - ox) - log(l - r) ' 



a 2 r 3 V J a? \{t + ox) 2 r 2 



1 W"')^ ' ^_f_U_ , 



a 2 r 2 V r b V 1-t J \(l + ^x) 2 
The equality (3.5) is equivalent to 



T — ^OO 



.2 



lim / e^ (x) = / e'^dx. (3.6) 
J(ft) 



By construction, we have 

<p T (0) = cp' T (0) = 0, c^(0) = -l, <Mx)<0, VieJ(S). 
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Let us observe that 



1 x 2 

lim ip T (x) = -<p"(0)x 2 = - — , \/x e M. (3.7) 



>o r ' y ' 2 rry ' 2 
Indeed, fix x € R and assume r is small enough so that 



t\x\ < -. (3.8) 



Observe that 



(i w 1 ^ 2(1 -r) <H ( a \ d? ' ( 1 



a 2 r 2 V r ^'"- u b V 1-r y dxi |X_U V (1 + fx) 2 
1 / 2(1 - r) / a Y , , , . / cr \ A 



1 -T 



Using the estimate cr = 0(r 2 ) as r -> we deduce that there exists C > such that, for any j > 
we have 

| #(0)| < C{j + l)Wi- 2 . 

Hence 

4lvi j) (ok j I < CjlrxP'-V, Vj > 2. 
Thus if t satisfies (3.8), we have 

Vr {x) + y = ^(s) - ^(0)* - i^(0)x 2 = £ i V W( y , 

i>3 J ' 

where the series in the right-hand side is absolutely convergent. Hence 

| ip T (x) + — | < Cx 2 \tx\ j\Tx\ j ~ 3 < C\Tx\x 2 ^2j2 j - 3 . 

j>3 j>3 

This proves (3.7). 

Next we want to prove that there exists a constant A > such that 

<Pt(x) < A(l - \x\), Vx£,J(h), Vr « 1. (3.9) 

We will achieve this by relying on the concavity of ip T over the interval J(h). The graph of ip T is 
situated below either of the lines tangent to the graph at x = ±1. Thus 

<p T (x) < p T (l) + y/ T (l)(x - 1) < -^ r (l) + ^ r (l)x, 

¥V(*) < Vr(-l) + ^(-1)(^ + 1) < ^(-1) + ^(-l). 

Now observe that 

d . , 1 / 2a 1 2cj 1 \ 2a ( 1 1 \ 

= -o-o -— ^ 5— + 



cix v ; a 2 r 2 I r 1-^ r (1 + f x) 3 i a 2 r 3 1(1 + fx) 3 l-^- x j' 

Using the fact that a = 0(t 2 ) we deduce from the above equality that 

|^ T (±1)|=0(1), asr^O. 
This proves (3.9). Using (3.7), (3.9) and the dominated convergence theorem we deduce 



T— >00 

We conclude that 



lim / e Vr ^dx = / e~^dx = V2tt. 



.,2 



p h ~ e* ft(T) a^\/2^ as h -> (3.10) 
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Now observe that 



, f v 1, h , 1 2(1 - r) log(l - r) 
4>h{r) = t log(l - r) 



/i ov ' a 2 r 2 a 2 r 2 t a 2 r 2 a 2 r 2 

Using (3.3) we deduce 



Also 



; , , 2(1— O 1 

3 </>ft(r) = (1 _ r )^3- e -^7 



In any case, using (3.3), (3.4) and (3.10) we deduce that 

logife(u>) ~ /clog a = Hog(c + 1) as k — > oo. (3.11) 

Thus 

log 5m (*>) = log = 7—r = 0, 

V Im+lM J 



so that 



lim q m = lim r m = 1. □ 

m— >oo m— >oo 



where 



Hence 



Example 3.5. If we let c = in the above example , then we deduce that 

r i i 

I k (w)= / t k e ^*dt ~ e MT) v / 2irv(h) 
Jo 

2 2 

Mr) ~ -3 [ " , v(h) ~ 1 (A X 3 

logI fc (u;) ~ -3 Q 

4. A PROBABILISTIC PROOF OF THE GAUSS-BONNET THEOREM 

Suppose that M is a smooth, compact, connected oriented manifold of even dimension m. For any 
Riemann metric g we can view the Riemann curvature tensor R g as a symmetric bundle morphism 
R g : A 2 TM ->■ A 2 TM. Equivalently, using the metric identification T*M = TM we can view /?„ 
as a section of A 2 T*M ® A 2 T*M. 

We will denote by QP'*(M) the sections of A P T*M ® A q T*M and we will refer to them of double 
forms of type (p, q). Thus R g G ri 2,2 (M). We have a natural product 

• : Q p ' q (M) x JF'' 9 '(M) -> r2 p+p ' ,9+9 '(M) 

defined in a natural way; see [1, Eq. (7.2.3)] for a precise definition. 

Using the metric g we can identify a double-form in J7 fc,fc (M) with a section of A k T*M® A k TM, 
i.e., with a bundle morphism A k TM — > A k TM and thus we have a linear map 

tr : n k ' k (M) -»• C°°(M). 

For 1 < < y we have a double form 

= R g » ■■■ • R g £ n 2k ' 2k {M). 
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We denote by dV g G m (M) the volume form on M denned by the metric g and the orientation on 
M. We set 

e 9 (M) := (27r) i (f) , tr(-^)dy g G O m (A/). 

The form e g (M) is called the Euler form of the metric g and the classical Gauss-Bonnet theorem 
states that 

/ e g (M) = x(M) =: the Euler characteristic of M. (4.1) 

In this section we will show that the Gauss-Bonnet theorem for any metric g is an immediate conse- 
quence of the Kac-Rice formula coupled with the approximation theorem Thm. 1.6. 

Fix a metric g. For simplicity we assume that vol 9 (M) = 1. This does not affect the generality 
since e cg (M) = e g (M) for any constant c > 0. Consider the random function u £ on M defined by 
(1.1, 1.2). Set 

.m+2 \ § 



Observe that for e > sufficiently small, any X, Y G Vect(M) and any p G M we have 

h £ (X(p), Y(p)) = E(Xv £ (p), Yv £ (p) ) 

where h £ is the metric on M that appears in the approximation theorem, Theorem 1.6. 

For any smooth function / : M — > K and any p G M we denote by HesSp(/) the Hessian of / at 
p defined in terms of the metric h £ . More precisely 

Hess^(/) = XYf(p) - (V £ x Y)f(p), \/X, Y G Vect(M), 

where V e denotes the Levi-Civita connection of the metric h £ . Using the metric h £ we can identify 
this Hessian with a symmetric linear operator 

Hess;(/) : (T p M,h £ ) -> (T p M,h £ ). 

For any p G M we have a random vector dv £ (p) G T*M. Its covariance form S(dv e (p)) is precisely 
the metric h e , and if we use the metric h e to identify this form with an operator we deduce that 
S(dv e (p)) is identified with the identity operator. 

For every smooth Morse function / on M and any integer < k < m we have a measure Vf^ on 

M 

d/(p)=0, ind(/,p)=fc 

where ind(/, p) denotes the Morse index of the critical point p of the Morse function /. We set 

m 
fc=0 

The Poincare-Hopf theorem implies that for any Morse function we have 

u f (dp)=x(M). (4.2) 

M 

Using the random Morse function v £ we obtain the random measures u v e p , u v e. We denote by vf. 
and respectively u £ expectations. The Kac-Rice formula implies that 

v k = ——mP%(p)\dV h e{p)\, 
(2vr)2 



L 
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where 



pUp) = det5 L(p) g (l det ^ s U v£ )\ I dv£ (p) = °> indHes 4(^) = fc 



= (-If El detHess^(v E ) | cfo £ (» = 0, indHess^ £ ) = kj. 

As shown in [1, Eq. (12. 2.11)], the Gaussian random variables Hess £ (i> £ ) and dv £ {p) are indepen- 
dent so that 



Thus 



pl(p) = (-l) fc £^detHess £ (?/) \ indHess £ (t> £ ) = k 
( 27r ) 2 t^o 

— ^£;(detHess £ (?; £ ) | indHess £ (u £ ) = k^j\dV h ,(p) 
E(detHess e p (v £ ))\dV h s(p)\. 



2 

k=0 



(27T) 

1 



(27T) 

Prom the Poincare-Hopf equality (4.2) we deduce 

X(M)= f u £ {dp) = —^f E(detHe SS £ p (v £ ))\dV h .(p)\. (4.3) 

Observe that Hessian Hess £ (/) of a function / can also be viewed as a double form 

Hess £ (/) G Q^M). 

In particular, Hess £ (?; £ ) is a random (1, l)-form and we have the following equality, [1, Lemma 
12.2.1] 

- 2R h e = ^(Hess £ (^ £ )' 2 ) , (4.4) 

where R^e denotes the Pdemann curvature tensor of the metric h £ . On the other hand we have the 
equality [ , Eq. (12.3.1)] 

detHess £ (u £ ) = — - trHess £ (^ £ ) ,m (4.5) 
ml 

Using (4.4), (4.5) and the algebraic identities in [1, Lemma 12.3.1] we conclude that 

1 / \ 1 / »^ 

—E[ detHessf, (v £ ) = m tr — R h i 

(2tt)t V pV >) (27r)f(f)! V h 

Using this equality in (4.3) we deduce 



X (M) = [ e h e{M), 

JM 



i.e., we have proved the Gauss-Bonnet theorem for the metric h £ . Now let e — > 0. As we have 
mentioned, Theorem 1.6 implies that h £ — > g so in the limit, the above equality reduced to the Gauss- 
Bonnet theorem for the original metric g. 
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Appendix A. Jets of the distance function 

Suppose that (M, g) is a smooth, m-dimensional manifold, p G M, U is an open, geodesically 
convex neighborhood of p and (x 1 , . . . , x m ) are normal coordinates on U centered at p . We have 
a smooth function 

77 : U X U ->• [0, oo), ^(p, q) = dist 9 (p, q) 2 . 
We want to investigate the partial derivatives of r at (p ,p ). Using the above normal coordinates we 
regard 77 as a function r\ = rj(x, y) denned in an open neighborhood of (0, 0) G M m x M m . 

If / = /(i 1 ) ■ ■ ■ ,t N ) is a smooth function defined in a neighborhood of G and is a 
nonnegative integer, then we denote by [/]& the degree fc-homogeneous part in the Taylor expansion 
of /atO, i.e., 

^ k = T\ E d?f\ t=0 t a eR[t\...,t N }. 

' \a\=k 

In the coordinates (x l ) the metric g has the form (using Einstein's summation convention throughout) 

g = gijdx l dx j , 
where satisfy the estimates [ , Cor. 9.8] 

9kt = hi ~ ^R ikj e(0Wx j + 0(\x\ 3 ). (A.l) 



We deduce that 



9 M = S ki + -i^(0)xV + 0(|x| 3 ). (A.2) 



1 

3 

The function rj satisfies a Hamilton- Jacobi equation, [27, p. 171], 

9 M d tf )d t/ =M*,y), Vx,y. (A.3) 
Moreover, 77 satisfies the obvious symmetry conditions 

m 

rj(x, y) = 7](y, x), 77(f), x) = r/(x, 0) = \x\ 2 := ^(x*) 2 . (A.4) 

i=l 

As shown in [7, Lemma 2.2] we have 

m 

lvh = \x- y \ 2 = Y,(* l -y 1 ) 2 - (A.5) 

i=l 

The symmetries (A.4) suggest the introduction of new coordinates (u, v) on U x [/, 

Uj = - y l , = a? 5 + y J . 

Then 

The equality (A.2) can be rewritten as 



J*( x ) =S M + ^J2 RikjM + Vi )( Uj + Vj ) + 0(3). (A.6) 
The symmetry relations (A.4) become 



12 



rj(u, v) = r](— u, v), i](u, u) = \u\ 2 , (A.7) 

while (A.5) changes to 

[77)1 = 0, [r]h = \u\ 2 . (A.8) 
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The equality (A. 3) can be rewritten 

9 k \x) ( r,' Uk + rf Vk ){ r,' ue +rf Vi )= 4 V - (A.9) 

k ' 1 =-A k =:A t 

Note that 

[A k } = [A e ] = [<? M ]i = 0, (A.10) 

while (A.8) implies that 

[A fc ]i = 2u k . 

We deduce 

4Ws = J>%( [A k ]i[Ai] 2 + [A fc ] 2 [B £ ]i) = J2 2 i A kh[A k ]i = 4^u k [A k ] 2 . 

k,£ k k 

We can rewrite this last equality as a differential equation for [77)3 namely 

[vh = ^2u k (d Uk +d Vk )[r]] 3 . 

k 

We set P = [77)3 so that P is a homogeneous polynomial of degree 3 in the variables u, v. Moreover, 
according to (A.7) the polynomial P is even in u and P(u, u) = 0. Thus P has the form 



P = Y J C i {u)v i +P (V), 



= :P 2 

where Cj(it) is a homogeneous polynomial of degree 2 in the variables u, and Pq(v) is homogeneous 
of degree 3 in the variables v. 
We have 

}^ u k d Vk P 2 = ^C fc (?j)nfc, Qi := y^u k d Vk Pp, ^7j fc 5 Ufe P = 0, 

Aj Aj 

and the classical Euler equations imply 
We deduce 

P = 2P 2 + Q 3 + Qi, 

where the polynomials Q3 and Qi are odd in the variable u. Since P is even in the variable u we 
deduce 

Qs + Qi = 0, 

so that P 2 + P = P = 2P 2 . Hence P 2 = P = and thus 

[77)3 = 0. (A. 11) 

In particular 

[Ajfc] 2 = 0, VA;. (A.12) 
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Going back to (A.9) and using (A. 10) and (A. 12) we deduce 

4W4 = Y^9 H h\Mi[Mi + YmilAkUAih + [A k ] 3 [Mi) 

(A.13) 

k,£ k 

We set P = [ry]4. The polynomial P is homogeneous of degree 4 in the variables u, v, and it is even 
in the variable it. We can write P = Pq + P 2 + P4, where 



Pa = ^ c ijk iUiUjU k U£, P 2 = Qij(u)viVj, 



and Po is homogeneous of degree 4 in the variables u, Qij(u) is a homogeneous quadratic polynomial 
in the variables ii. We have 

fc fc 

We have 

^ u k d Uk P 2v = 2uP 2u , v = 0, 1, 2. 
^u fc <% fc P 4 = 0, 

y]u fc 9„ fc P 2 = ^UkQijidkiVj + 5 k jVi) = ^(QkjUkVj + Qj k VjU k ) 

k k,i,j k,j 

Using these equalities in (A.13) we deduce 

4P 4 + 4P 2 + 4P = 4 Y,\9 k %u k u t + 4P 4 + 2P 2 + ^ u k d Vk P 

k,£ k 

+ ^2[Qjk + Qkj)ukVj- 

This implies Pq = so that P = P4 + P 2 , and we can then rewrite the above equality as 

P 2 = 2 ^[g M ] 2 u k u t + ^2{Q jk + Qkj)ukVj. (A.14) 

k,l k,j 

Note that the equality r(u, u) = |u| 2 implies P(u, u) = so that 

Pi(m) = Pt(«,w) = -P 2 (u,u). 

There fore it suffices to determine P 2 . This can be achieved using the equality (A.6) in (A.14). We 
have 

^^[g^hUkUi = - ^ Ri k ji(Ui + Vi)(Uj + Vj)u k Ul 

= g (X! R ikj£U k U<> ) + ^ Sj(u)vj, 

i,j k,l j 
* — v ' 
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where Sj(u) denotes a homogeneous polynomial of degree 3 in u. The equality (A. 14) can now be 
rewritten as 

^ Qij(u)viVj = - ^2 Qij{u)viVj + ^2 Sj(u)vj + - ^(Qjk + Qkj^jukVj- 

i,j i,3 3 k,j 

From this we read easily 



This determines Po. 



P2{u,v) = -^2 l Q i j(u)v i v j . (A.15) 



6 



As we have indicated above P 2 determines P4. 

P A (u) = -P 2 {u,u) = -- ^2 Rikj£UiUjU k u e . (A.16) 



6 

i,j,k,£ 

The skew symmetries of the Riemann tensor imply that P4 = so that 

[r]} A (u,v) = ^^2Qij{u)viVj, Qij(u) = y^RjkjtUkUz. (A.17) 

i,j k,l 

Example A.l. Suppose that M is a surface, i.e., m = 2. Set 

K = i?1212 = -^2121 = — -^1221- 

Note that K is the Gaussian curvature of the surface. Then 

Qll = ^2 R lkltUkUt = Kul, Q22 = ^2 R 1k2tUkUt = Ku\. 



)\2 = ^2 Rlk2lUkUl = ~ Ku l u 2 = Q21- 



k.C 



Hence 



P2(u,v) = ^r{u 2 2 vl + u\v\ - 2uiu 2 viV2) = ^-(uiv 2 - u 2 vi) 2 . 
o o 



□ 



Appendix B. Spectral estimates 
As we have already mentioned, the correlation function 

<? £ (P,l) = Y, w ^^^k(p)^k(q) 

k>0 

is the Schwartz kernel of the smoothing operator w £ (\^A). In this appendix we present in some detail 
information about the behavior along the diagonal of this kernel as e — > 0. We will achieve this by 
relying on the wave kernel technique pioneered by L. Hormander, [18]. 

The fact that such asymptotics exist and can be obtained in this fashion is well known to experts; 
see e.g [i i] or [29, Chap.XII]. However, we could not find any reference describing these asymptotics 
with the level of specificity needed for the considerations in this paper. 
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Theorem B.l. Suppose that w G S(R) is an even, nonnegative Schwartz function, and (M,g) is a 
smooth, compact, connected m-dimensional Riemann manifold. We define 

<T :MxM^R, £ e {p,q) = Y,™(^k)^k(p)*k{q), 



k>0 



where (^*k)k>i is an orthonormal basis of L (M,g) consisting of eigenf unctions of A g . 

Fix a point p G M and normal coordinates at p defined in an open neighborhood Oq of p . 
The restriction of <§ e to <§ e to Oo x Oo can be viewed as a function t§ e (x,y) defined in an open 
neighborhood of (0,0) inR m x R m . Fix multi-indices a,/3 G (Z>o) m . Then 



d^%x,y)\ x=y=0 = s- m - 2d ^ t -^ n^w(\x\)x a+ ?dx + 0(e 2 A , e -+ 0, (B.l) 
where 

d{a,P):= — - — 

Moreover, the constant implied by the symbol 0(e) in (B.l) uniformly bounded with respect to p . 

Proof. For the reader's convenience and for later use, we decided to go in some detail through the 
process of obtaining these asymptotics. We skip many analytical steps that are well covered in [20, 
Chap. 17] or [26]. 

Observe that for any smooth f : M ->lwe have 

w £ (VA)f = ±- f w £ (t)e u ^fdt = -L / w (-) e^fdt. (B.2) 

The Fourier transform w(t) is a Schwartz function so w(t/e) is really small for t outside a small 
interval around and e sufficiently small. Thus a good understanding of the kernel of e ltv ^ for t 
sufficiently small could potentially lead to a good understanding of the Schwartz kernel of w £ (\^K). 
Fortunately, good short time asymptotics for the wave kernel are available. We will describe one such 
method going back to Hadamard, [17, 27]. Our presentation follows closely [20, §14.4] but we also 
refer to [26] where we have substantially expanded the often dense presentation in [20]. 

To describe these asymptotics we need to introduce some important families homogeneous gen- 
eralized functions (or distributions) on R. We will denote by C~°°(r2) the space of generalized 
functions on smooth manifold Q,, defined as the dual of the space compactly supported 1-densities, 
[15, Chap. VI]. 

For any a G C, Re a > 1 we define x+ : M — )■ R by 

X a +(x) = r , 1 ^ 4' x + = max(x,0). 
T(a + 1) 

Observe that we have the following equality in the sense of distributions 



^xV 1 = X a +(x), Rea>l. 



We can use this to define for any a G 



X+ := ^Xt k e C-°°(R), fc>l-Rea. 



For Re a > we denote by \x\ a the generalized function defined by the locally integrable function 

1 



lxl°0»0 = t^+ttW 



2 
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The correspondence a i— > |x| a is a holomorphic map {Re z > 0} — > C~°°(K) which admits a holo- 
morphic extension to the whole complex plane, [14, Chap. 1], [26]. This is a temperate generalized 
function, and its Fourier transform is given by, [14, 26], 

\xT a (6 = ^2 a+1 | X r (a+1) (0, Va G C. (B.3) 

Denote by K t (x, y) the Schwartz kernel of e lt ^. We then have the following result [20, §17.4] 
or [26]. 

Theorem B.2. Set n := m + 1, and let 

r](x, y) = dist fl (x, yf, x, y G M. 

There exists a positive constant c > 0, smaller than the injectivity radius of (M,g), such that for 
dist 9 (x, y) < ewe have the following asymptotic expansion as t — >■ 

oo 

K t (p,q)~J2 U k(P,Q)d m (2k)M k {t,p,q), \t\ < c, (B.4) 

k=l 

where for Re a > we have 

X + 2 {tl-r)( P) q))- X+ 2 {t 2 _- V {p,q)) 

r(2^i) 



d m (2a) 



vr^r(2a) 



Let us explain in more detail the meaning of the above result. The functions Uk are smooth func- 
tions defined in the neighborhood dist s (p, q) < c of the diagonal in M x M. For fixed q, the 
functions p \- > Vk{p) '■= Uk{p, q) are determined as follows. 

Fix normal coordinates x at q, set \g\ := det(<7y), and 

h( x ) : = -^(Vlogls-I, x) = -^^g^x^d^ ]og\g\. 

Then Vk(x) are the unique solutions of the differential recurrences 

Vi(0) = 1, 2x ■ VVi = hV 1: \x\ < c, (B.5) 

^x-VV k+1 +(^l-^}?jV k+1 = -A g V k , V k+ i(0)=0, \x\<c, k>l. (B.6) 

We have the following important equality 

lim X a (t,p,q) = \x\ 2a ~ 2 ~ m (t), VaeC. (B.7) 

dist s (p,<jf)— >0 

The asymptotic estimate (B.4) signifies that for any positive integer n there exists a positive integer 
N(p) so that for any N > N(fi) the tail 

TV 

7 N {t,p,q) := K t (p,q) -^U k (p,q)d m {2k)'K k (t,p,q) 

k=l 

belongs to C At ( (— c, c) x M x M ) and satisfies the estimates 

||^T^(i,-,-)|| CM _ J(MxM) < Cltl 2 "-"- 1 -", |t|<c, ,j<fi, N>N(ji). (B.8) 

Fix a point p G M and normal coordinates at p defined in a neighborhood Oo of p . Then we 
can identify a point (p, q) G Oo x Oo with a point (a;, y) in a neighborhood of (0, 0) in 
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Using (B.2) we deduce 



(B.9) 



:K 



a,0 



Choose an even, nonnegative cutoff function p G (R) such that 



tit) 



i. 1*1 <i 

v 0, |t|>f, 
where c > is the constant in Theorem B.2. Then 



d a x d^{x,y)\ x=y = - £ (K^,p(t)w (?) ) + - £ (k^, ( 1 - p(t))i 
Let us observe that that for any N > 

±(K^,(l-p(t))w(^)=0(e N ) ase^O 

ViV > d-d^%x,y)\ x=y ~ i(^,p(t)^ 0) ) +0(e JV ), e ->• 0. 



Thus 



On the other hand 



Recall that 



d^K t (x,y) ^J2d m (2k)d^{U k (x,y)^ k (t,x,y)} 



(B.10) 



(B.ll) 



k=l 



d{a,P) 



\a + P\ 



One can show (see [7, 26]) 



<%9yKt(x,y)\x=y=o ~ ^Z" 4 ™'^ 



ftfclXl 



-m-2d(a,l3)+2k 



(B.12) 



fc=0 



where A m)Q is a universal constant depending o«(y on m, a, /?, which is equal to if \a + /3\ is 
odd. 



Lemma B.3. (a) For any r £ Z owe? any N > Owe have 

±(\x\ r ,pw £ ) = e r ((\ X \ r ,w) + 0(e N )) 
(b) For every positive integer r we have 

/I I— t ^\ V / i \r — 1 

(Ixl >«>)- 



as e — > 0. 



r(5) 



t\ t w{t)(1t. 



Proof, (a) For transparency we will use the integral notation for the pairing between a generalized 
function and a test function. We have 

{\x\ r ,We) = - ! \x\ r (t)p(t)w(t/e)dt= [ \ x \ r (et)p(et)w(t)dt 



\ X \ r (t)p(et)w(t)dt = e r {\x\ r ,p e w), pS) = 
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Now observe that p £ w — w = w(p £ — 1) — )• in S(R). More precisely for k > we have 



(p £ - l) = Q{e N t N ) as e -> 



dk 



This implies that 
so that 
(b) We have 



(\ X \ r ,w(p £ -l))=0(e N ) ase^O, 

(lx| r , Pew) = (\x\ r , w)+{\ X \ r , - 1) > = (lx| r , w> + 0{e N ) as 0. 

(B.3) 

r| r_i i/;(T)dT. 



(\xr,w) = (bP r ^) = v^-'iixr^Mr)) 

_ f u,r-l„ 



□ 

Using (B.10) and the above lemma we deduce 

(B.13) 

where D m ^p is a universal constant that depends only on m, a, (3 which is = if |a + /3| is odd, 



r(|) 



r| r_1 'w(r)dr, r = m + 2d(a, /5). 



(B.14) 



To determine the constant D m a ^ it suffices to compute it for one particular m-dimensional Rie- 
mann manifold. Assume that (M, g) is the torus T m equipped with the flat metric 

m 

g = J2(d6 1 ) 2 , < 9 i < 2tt. 
i=l 

The eigenvalues of the corresponding Laplacian A m are 

\k\ 2 , k = (h,...,k m )ez m . 

Denote by -< the lexicographic order on 7L m . For 6 = (8 1 , . . . , 9 m ) £ K and fc e Z m we set 

m 

A real orthornormal basis of L 2 (T m ) is given by the functions 



Then 



so that 



(27T) 



1, fc = 

25 S in(M), fc^O, 
cos(k,6), k<0. 



1 



(27T) 



9-0) 



fceZ" 
;H-|/3| 
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Define 

W m ,u £ :R m ^R, W m (x)=w(\x\), u £ (x) = W m (ex)x a+IB . 
Using the Poisson summation formula [19, §7.2] we deduce 

Observe that 

u £ (0= [ e- i ^ x) w{e\x\)x a+ Pdx = {id i ) a+ P ( [ e~ i{ ^ x) W m {ex)dx 



Hence 

As e — > we have 

02^(0,0) = ^ m " |a+/3 '^^ ( (*%) Q+ ^Wm(0) + 0(e*)) , VJV. 
Now observe that 

(ia 5 )" +/3 En(0) = / W(\x\)x a+I3 dx. 

so that 

95^(0,0) = £ - m -l«+/3|^L^. ^^IK+^ + O^)) , VAT. 
This shows that 

{ Z7r J JR m 

This completes the proof of Theorem B.l. 

Remark B.4. Note that 

/ w(\x\)x a+ Pdx = \ [ x a+P dA{x) \ ( r w{r)r m+ \ a+ ^- l dr 
jR m \J\x\=1 J \Jo 

On the other hand, according to [ , Lemma 9.3.10] we have 



/ x a+ ?dA{x) = Z m ^ := < 

J\x\=\ 



'1*1=1 

We can now rewrite (B.16) as 



( 9rr fe r( c *i + Pi +1 '\ 
0, otherwise. 



L>m,<x,P=£ ' WW' 
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Proposition B.5. Fix a point p £ M and normal coordinates {x l ) near p. For i ^ j we denote by 
K{j (p) the sectional curvature of g at p along the plane spanned by d x i , d x j . For any multi-induces 
a,/3£ (Z> ) m weset 

^ :=d^%x,y)\ x=y=Q . 
Then there exists a universal constant Z m that depends only on the dimension of M and the weight 
w such that 

g kii ~ S tm = Z ™ K *J (P) £ ~ m ~ 2 (! + °( £2 ) ) as e^Q. (B.19) 



Proof. Using (B.12) we deduce 

~ ~ \{K?" - Kr\ v (t)w (~) ) + 0(e N ), e -+ (B.20) 
On the other hand from (B.9) we conclue 

oo 

- ~ d ^ 2k ) ~ ^ rKn) { U ^ yWkit, x, y) } \ x=y=0 (B.21) 

k=l 

To investigate the above asymptotics we use the technology in [26]. 

Let us introduce some notations. For a positive integer k we denote by d k a generic mixed-partial 
derivative of order k in the variables x % , yK We denote by d k 7] the collection of A;-th order derivatives 
of r]{x,y). 'Pi(X) will denote a homogeneous polynomial of degree i in the variables X, while 
J > / c (X)J > £(y) will denote a polynomial which is homogeneous of degree k in the variables X and of 
degree I in the variables Y. We then have the equalities 

■K a = ViidrfiHa-u (B.22) 
8 2l K a = y 2 {di 1 )'K a _ 2 + , ? 1 {d 2 r ] )'K a - 1 , (B.23) 
d*K a = ? 3 (dr,)H a - 3 + 5 i(ar ? )T 1 (a 2 r ? )?f a _ 2 + JiC^JCa-i, (B.24) 

o 4 K a = y i {di 1 )'K a _ A + (T 2 (a7?)yi(a 2 7?))?{ a _3 

To simplify the presentation we will assume that 

i = l, j = 2. 

Also we will denote by O(l) a function f(x,y) such that f(x,y)\ x=y =o = 0. The computations in 
Section A show that for x = terms of the form < ?j{drj) and < J > k (d 3 ?]) are O(l). In particular, the 
above equalities show that the 1st and 3rd order derivatives of "K a are O(l). We have 

3%iefr(u k M k ) = a 2 ! ( {d 2 y2 u k )x k + 2d y 2U k d y2 j{ k + u k d^ k ) 

= {d 2 xl d 2 y2 U k )X k + (fl£l7 fc ) + {dfrUk) (d 2 2 M k ) (B-26) 

+4(d 2 v c4)(3 2 v iK fc ) + u k d 2 xl d 2 y2 -K k + 0(1), 
<9 2 lx2 d 2 v (C4iK fc ) = d 2 xlx2 ( (d 2 v c4)?t fc + d y iU k d y2 K k + d y 2U k d y iX k + ^a 2 lj/2 ?c fc ) 

= (^ 2 ^ ly2 C/ fe )^ + (^^(^Wfc) 

+dl 2y iU k dl ly2 Ji k + d 2 ^U k d x \ y% 'K k + dl 2y2 U k d x i y i^{ k + d x \ y2 U k d x \ yX 'K k 

+d 2 x i x2 U k d 2 yly2 X k + ^<9 2 lx2 d 2 v Jt fc + 0(1) 

(B.27) 
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He deduce that 

4 

(d 2 xl o 2 y2 - d 2 xlx2 d 2 yly2 )(u k x k ) x=y=0 = Y,nx k -i\ x=y=0 , 

j=Q 

where the coefficients Tl are polynomials in the derivatives of U k and r\ at (x,y) = (0,0). Using 
(B.22)-(B.25) we deduce 

n = n = o. 

Moreover, the terms that appear in T? appear only when we take forth order derivatives of "K k . 
Upon inspecting (B.26) and (B.27) we see that the 4th order derivatives of < K k are multiplied by U k . 
According to (B.6) the function U k is O(l) if k > 1. Hence = for k > 1. We deduce 

oo 

K ii,jj _ K m „ J- d m (2k) {T°X k + T k "K k _i + T 2 X k _ 2 ) \ x=y=0 

k=l 

= B-x3i-i\ x= y = o + Bo9io\ x=y= o + BiKi\ x=y= Q + • • • , 

where 

B_i = d m {2)Tf, B = d m {2)Tl Bt = d m (2)T° + d m (4)T}, . . . . 
The term B-\ can be alternatively described as 

B—\ — A- m ii-jjQ ^4m,ij';jj,0) 

where the coefficients A ma a are defined as in (B.12). Using (B.14) and (B.16) we deduce 

5_i = 0. 

To compute If we observe first that 

7](x — y) = '^2(x i — y 1 ) 2 + higer order terms. (B.28) 

i 

Using (B.23) we can simplify (B.26) and (B.27) in the case k = 1 as follows. 

^^TO) = (^xfigat^fti + Ujd^d 2 ^ + 0(1), (B.29) 



(B.30) 



+^ aya cr l a B i v iiKi + t^i^v^ 1 + o(i). 

Using (B.23), (B.25) and (B.28) we deduce that 

Tl = (d 2 xl d 2 2 -d 2 xlx2 d 2 yly2 )r,\ m 
+2 (d^ + ^ 3 l7i) | ( o,o) + 2 (^i v iCTa + d^U,) \ m . 

Using the transport equation (B.5) we obtain as in [10, VI.3] that U\ coincides with the function 
ip(x,y) in [10, VI.3 Eq.(33)] or the function uo(x,y) in [6, p. 380]. For our purposes an explicit 
description of U\ is not needed. All we care is that 

Ui(x,y) = U 1 (y,x), Ui{x,x) = l. 

These conditions imply that the Hessian of U\{x,y) at (0,0) is a quadratic form in the variables 

U{ = {x l — y l ) so that 

^ 1 C7 1 (0,0) + ^ v C/ 1 (0,0) = ^ 2 C7 1 (0,0)+^ v i7i(0,0) = 0. 

Hence 

n = (a 2 xl d 2 y2 - e^v) iko,o)- 
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Using (A. 17) we conclude that 

Tl = ZR 1212 = ZK 12 (p), 
where Z is a universal constant, independent of (M, g). Hence 

K?' jj - Kl j ' lj ~ d m (2)ZK 12 (p)X°\ x=y=0 + Y, B ^k\ x=y =o- 

The equality (B.19) now follows from the above equality by using (B.20), (B.7) and Lemma B.3. □ 

Appendix C. Gaussian measures and Gaussian vectors 

For the reader's convenience we survey here a few basic facts about Gaussian measures. For more 
details we refer to [8]. A Gaussian measure on R is a Borel measure 7^ )W , v > 0, m G R, of the form 

i \ 1 -(— "? a , 
T^W = ^ — e a« dx. 

The scalar /z is called the mean, while u is called the variance. We allow v to be zero in which case 

7/i,o = Sfj, = the Dirac measure on R concentrated at p. 
For a real valued random variable X we write 

XeN(fi,v) (CI) 

if the probability measure of X is 7^. 

Suppose that V is a finite dimensional vector space. A Gaussian measure on V is a Borel measure 
7 on V such that, for any £ 6 V v , the pushforward (7) is a Gaussian measure on R, 

£*(7) = Mt)At) ■ 

One can show that the map V v 9 £ 1— >• //(£) G R is hnear, and thus can be identified with a vector 
/x 7 G V called the barycenter or expectation of 7 that can be alternatively defined by the equality 



H = J vdj(v). 
Jv 



Moreover, there exists a nonnegative definite, symmetric bilinear map 

S : V y x V v -> R such that = S(£,£), V£ G V y . 

The form S is called the covariance form and can be identified with a linear operator S : V y — > V 
such that 

E(£,i7) = (£,Sfy>, V£, ?? GF V , 

where (— , — ) : V v x V — > R denotes the natural bilinear pairing between a vector space and its 
dual. The operator S is called the covariance operator and it is explicitly described by the integral 
formula 

{£,Sri) = £(£,r?) = y - »y)(v,v- /x 7 > 

The Gaussian measure is said to be nondegenerate if 5] is nondegenerate, and it is called centered if 
/x = 0. A nondegenerate Gaussian measure on V is uniquely determined by its covariance form and 
its barycenter. 
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Example C.l. Suppose that U is an n-dimensional Euclidean space with inner product (— , — ). We 
use the inner product to identify U with its dual U y . If A :U—> U is a symmetric, positive definite 
operator, then 

~f A (du) = e-^A-^u) | d , (C 2) 

' AK ' (2vr)fVd^4 11 
is a centered Gaussian measure on U with covariance form described by the operator A. □ 

If V is a finite dimensional vector space equipped with a Gaussian measure 7 and L : V — > U is 
a linear map, then the pushforward L^j is a Gaussian measure on U with barycenter 

and covariance form 

£l» 7 : U y x £/ v -> K, S L , 7 (? ? ,r ? ) = S 7 (Z v r/,L v r ? ), V?] € C/ v , 

where L v : £/ v — > V v is the dual (transpose) of the linear map L. Observe that if 7 is nondegenerate 
and L is surjective, then ^7 is also nondegenerate. 

Suppose (S, n) is a probability space. A Gaussian random vector on (S, /u) is a (Borel) measurable 
map 

X : S — > V, V finite dimensional vector space 

such that ,X*/z is a Gaussian measure on V. We will refer to this measure as the associated Gauss- 
ian measure, we denote it by "fx an d we denote by (respectively S(X)) its covariance form 
(respectively operator), 

Vx(Zi,&) = E({Z 1 ,X-E(X)){£ 2 ,X-E(X))). 

Note that the expectation of jx is precisely the expectation of X. The random vector is called 
nondegenerate, respectively centered, if the Gaussian measure 7^ is such. 

Let us point out that if X : S — > U is a Gaussian random vector and L : U — > V is a linear map, 
then the random vector LX : S — > V is also Gaussian. Moreover 

E(LX) = LE(X), V LX (Z,0 = V x (L v !i,L v Z), V£ 6 V v , 

where Z v : F v -»• U v is the linear map dual to L. Equivalently, S(LX) = LS(X)L V . 

Suppose that Xj : S — > V\, j = 1, 2, are two centered Gaussian random vectors such that the 
direct sum X\ © X2 : S — > V\ © V2 is also a centered Gaussian random vector with associated 
Gaussian measure 

lXx®x 2 = PXi®x 2 ( x i, x 2 )\dx\dx 2 1. 

We obtain a bilinear form 

cov(X 1 ,X 2 ) -.VXxV^R, cov(X 1 ,X 2 )(ti,&) = £(6,6), 

called the covariance form. The random vectors X\ and X2 are independent if and only if they are 
uncorrelated, i.e., 

cov(Xx,X 2 ) = 0. 

We can then identify cov (X±, X 2 ) with a linear operator Cov(X\, X 2 ) : V 2 — > V\, via the equality 

E((^ 1 ,X 1 )(^,X 2 )) = coviX^X^,^) 

= (Z 1 ,Cov(x 1> x 2 )&), vaeVi 7 , e 2 e^, 

where £} 2 G V2 denotes the vector metric dual to £2- The operator Cov(X\, X 2 ) is called the 
covariance operator of X\, X 2 . 
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The conditional random variable (Xi\X2 = £2) has probability density 

PXiex 2 (a5i)«2) 



P(X 1 \X2=x 2 ){ x l) = 7 7 \Tj 1 

v 1 ' J Vl Px 1 ®x 2 (xi,X2)\dx 1 \ 



For a measurable function / : V± — > M the conditional expectation £(/(Xi)|X2 = X2) is the 
(deterministic) scalar 

£(/(Xi)|X 2 = £c 2 ) = / /(aJiWjXa^CaJiJIdxil. 

If X 2 is nondegenerate, the regression formula, [5], implies that the random vector (Jfil^ = x 2 ) is 
a Gausian vector with covariance operator 

S(Y) = SpGj - Cov(X 1 ,X 2 )S(X 2 y 1 Cov(X 2 ,X 1 ), (C.3) 

and mean 

E(X 1 \X 2 =x 2 ) = Cx 2 , (C.4) 

where C is given by 

C = Cov(X 1 ,X 2 )S(X 2 r 1 . (C.5) 

Appendix D. A class of random symmetric matrices 

We denote by Sym m the space of real symmetric mxm matrices. This is an Euclidean space with 
respect to the inner product 

(A,B) :=tr(AB). 

This inner product is invariant with respect to the action of SO(m) on Sym m . We set 

ft _\ E iv i = j 

VV2 i ^ <J - 

The collection (Eij)i<j is a basis of Sym m orthonormal with respect to the above inner product. We 
set 

\\f2aij, i<j. 

The collection (a,ij)i<j the orthonormal basis of Sym^ dual to (Eij). The volume density induced 
by this metric is 

\dA\ := YldSij = 2s(a) JJ day. 

Throughout the paper we encountered a 2-parameter family of Gaussian probability measures on 
Sym, m . More precisely for any real numbers u, v such that 

v > 0, mu + 2v > 0, 

we denote by Sym^ the space Sym m equipped with the centered Gaussian measure dF u>v (A) 
uniquely determined by the covariance equalities 

E(aijake) = u5ij5ke + v(5ik5je + S^Sjk), VI < i,j, .k,£ < m. 

In particular we have 

E(a% i ) = u + 2v, E(auajj)=u, E(d% j )=v, VI < i / j < m, 

while all other covariances are trivial. The ensemble Sym^ is a rescaled version of of the Gaussian 
Orthogonal Ensemble (GOE) and we will refer to it as GOE^. 
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For u > the ensemble Syrn^ can be given an alternate description. More precisely a random 
A G Sym^ u can be described as a sum 

A = B + Xl m , B G GOE^, X G N(0, u), B and X independent. 

We write this 

Symr/ = GOE^+AT(0,n)l m , (D.l) 
where + indicates a sum of independent variables. 

The Gaussian measure dT U)V coincides with the Gaussian measure dT u+ 2v,u,v defined in [ 14, App. 
B]. We recall a few facts from [ , App. B]. 

The probability density dT UjV has the explicit description 

where 



m(m+I) 

(2tt)^^ v / ^(^) 



D(u, v) = (2v) {m - 1)+ ^) (mu + 2v), 

and 



/ If 1 1 

it = — 



m \mu + 2v 2v J 2v(mu + 2v) 
In the special case GOE^ we have u = u' = and 

dToAA) = \^e-h**\ d A\. (D.2) 

(2irv) 4 

We have a Weyl integration formula [3] which states that if / : Sym m — > Ris a measurable function 
which is invariant under conjugation, then the the value f(A) at A G Sym m depends only on the 
eigenvalues Ai(^4) < ■ ■ ■ < X n (A) of A and we have 



E GO w m (f(x))=-± 7 - f /(Ai, . . . ,x m ) ( n i^-A^ijn 



_>4 

e 4« \d\\ ■ ■ ■ dX r 



— '-Qm,v (A) 

where the normalization constant Z m (v) is defined by 



(D.3) 



Z m (v)= / |Aj - Xj\ Y[e~^\d\i ■ ■ -dX m \ 

/TOm 

,yM l<«<j<m 2=1 

/m 2 
1 — r i — r *\' 
|A i -A i |TTe-^|dAi---dA r 



l<i<j'<m i=l 



The precise value of Z m can be computed via Selberg integrals, [3, Eq. (2.5.1 1)], and we have 

m 



Z rn = (2vr)f ml f[ = 2 f m! (f ) . (D.4) 



For any positive integer n we define the normalized 1-point corelation function p n:V (x) of GOE^ to 
be 

Pn,v{x) = r . . . / Qn^^j A 2 , • • • , A n )dAi • • • dX n . 

ZnVP) J R n-l 
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For any Borel measurable function / : R — > R we have [11, §4.4] 

-E GOKi {trf(X))= [ f(X)p n>v (X)d\. (D.5) 
The equality (D.5) characterizes p n>v . Let us observe that for any constant c > 0, if 

A G GOE^ ^cA G GOE^ v . 
Hence for any Borel set B C R we have 



L 



Pn^v{x)dx = I p thV (y)dy. 

IcB JB 

We conclude that 

CPn,c^v(. c y) = Pn,v{y)i Vn ' C ' V- ( D - 6 ) 

The behavior of the 1-point correlation function p n ,v(%) f° r n large is described by Wigner semicircle 
law which states that for any v > the sequence of measure on R 

p nvn -i(x)dx = n.2 p ntV (n^x)dx 
converges weakly as n — > oo to the semicircle distribution 

poo jV (x)\dx\ = I { \ x \< 2 ^ } ^\/ Av - x 2 \dx\. 

The expected value of the absolute value of the determinant of of a random A G GOE^ can be 
expressed neatly in terms of the correlation function p m+ i „. More precisely, we have the following 
result first observed by Y.V. Fyodorov [13] in a context related to ours. 

Lemma D.l. Suppose v > 0. Then for any c G R we have 

/ , , . \ i \ a , m+i / m + 3 \ c 2 

^GOE«,(|det(A-cl m )|) = 22(2?;) — T — — e^p m+1:V (c). 



Proof. Using the Weyl integration formula we deduce 

1 f m A 2 

^GOE^(|det(A-cl m )|) = - / Te-*|c-A i |TT|A i -A i |dAi---dA„ 

Z m {v) Jwn , =1 £A 

/ TT e ~4^ |c - Aj| 1 [ |Aj - A,|dAi • • • d\ m 



Z m (u) Z m+1 (y 



Qm+i.ufe Ai, . . . , A m )dAi • • • d\ Ti 



c c 

e- iv Z m+ l(v) rn + ieiv Z m+ i 

-Pm+l,v{C)=V 2 Pm+l,v(C) 



/ x /— / "i+i c 2 / m + 1 \ , . 3, , m+i _ / m + 3 A , . 

(m + 1)V2(2«)T- e ^r — — p m+ i,„(c) = 25(2?;)— r — — e5 %+1 ( c ) 



□ 



The above result admits the following generalization, [4, Lemma 3.2.3]. 
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Lemma D.2. Let u,v > 0. Set 

9 ti,v( x ) ■= Pm+l,v(x)e 

Then 



„ / , , / , , . \ 3 , , m + 1 _ f m + 3\ 1 /" , . (c~x) 2 x 2 , 



(D.7a) 

2 f(2t,)^rfet* > ) (7«*^ +1 ,J(c). (D.7b) 



In particular, ifu = 2kv, k < 1 we have 



o /m-L^X 1 /• 1 fa I f 2 c) 2 I (t fe +1)c2 

S Svm2fcu , u (|det(A- C l m )|) =22( 2w )fr (— ±i -= / p m+M ( C -x)e ^ fc 4 " cfe, 



(X := c — x) 



where 



2*(2v)TT[—j-\-j= I p m+ i,,(A)e 



Proof. Recall the equality (D.l) 

SynC" = GOE^+AT(0,u)l m . 

We deduce that 

^s y m- ( I det(,4 - cl m )| ) = S(det(S + (X - c)l)| ) 



= -j=jE GO EZ t {\det(B-(c-X)l m )\ \ X = x)e 

If 2 
= -= / £ GOE ^(|det(£ -{c-x)t m )\)e-^dx 

2/ , m+l _ / m + 3\ 1 f . . (c-j:) 2 x 2 

22(2u) 2 rl — - — j-j==J p m+1>v (c-x)e 4, 2uC fa;. 



Now observe that if u = 2/ct; then 



(c X X 1 . o o. 

^ = — H (x 2 - 2cx + c 2 ) 

4u 2u Akv Av 



□ 
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